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(Dated: November 5, 2009) 

The effect of the strong intersite Coulomb correlations on the formation of the electron structure of 
the t — V-model has been studied. A qualitatively new result has been obtained which consists in the 
occurrence of a split-off band of the Fermi states. The spectral intensity of this band increases with 
the enhancement of a doping level and is determined by the mean-square fluctuation of the occupation 
numbers. This leads to the qualitative change in the structure of the electron density of states. 
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The key point of the theory of the strongly correlated 
electron systems is the statement on a principal role of the 
one-site Coulomb repulsion of two electrons with the op- 
posite spin projections, e.g., the Hubbard correlations 111], 
in the formation of the ground state and the elementary ex- 
citation spectrum Q]. One of the brightest manifestations 
of the Hubbard correlations is splitting of the initial band 
of the energy spectrum into two Hubbard subbands when 
one-site energy of the Coulomb interaction U exceeds the 
bandwidth W. 

As was noted previously Jl 0], when hole concen- 
tration in a system described by the Hubbard model is 
small and value U is large the intersite Coulomb interac- 
tion starts playing an important role due to its relatively 
weak screening. In this case at distances close to inter- 
atomic the order of magnitude of characteristic value V 
of the Coulomb interaction between electrons can be com- 
parable with that of value U. Manifestation of the in- 
tersite correlations in the physical properties of systems 
has been considered in many works (see, for instance, 

II S ll & B S II El El El El El El). However, 
until recent time no attention has been paid to the fact that 
the strong intersite correlations (SIC) can cause splitting 
energy bands of the Fermi states into subbands 111 611 . 

Qualitatively, the physical origin of this phenomenon is 
similar to that of the occurrence of the Hubbard subbands 
due to the strong one-site correlations and is related to the 
fact that when the intersite Coulomb interactions are taken 
into account the energy of the electron located on site / be- 
comes dependent of the valence configuration of its nearest 
neighborhood. For example in the case of a lower Hubbard 
subband, if there is one electron on each of the z nearest 
neighboring ions, the "setting" energy of an electron is de- 
termined as a sum of one- site energy e and energy zV of 
the interaction. If a hole appears in the nearest neighbor- 
hood, the "setting" energy is e + (z — 1)V, i.e., its value 
is less than the previous one by V. 

These simple arguments show that when configuration 
neighborhood deviates from nominal one, one should ex- 
pect the occurrence of the states with the energies different 
by V in the Fermi excitation spectrum. For the strongly 
correlated systems with a relatively narrow energy band 



of the Fermi states, the situation becomes possible when 
hopping parameter t is commensurable with or less than 
V. Under these conditions splitting the initial band of the 
Fermi states is expected. Obviously, the more probable the 
deviation of the electron configuration of the neighborhood 
from the nominal one the higher is the spectral intensity of 
the split-off band. A quantitative measure of this deviation 
is a mean-square fluctuation of the occupation numbers. 
For this reason the split-off band was named the band of the 
fluctuation states (BFS) L16]. For undoped Mott-Hubbard 
insulators the electron configuration of the neighborhood 
corresponds to the nominal one. Upon doping the deviation 
occurs and in the energy spectrum of the Fermi states the 
BFS appears whose spectral intensity increases with dop- 
ing level. The aim of this Letter is to prove the presented 
qualitative interpretation. 

In order to demonstrate the effect as brightly as possible 
we will consider the Hubbard model in the regime of the 
strong electron correlations (U = oo) at electron concen- 
trations n <C 1. In this case the electron properties of the 
model are determined by the lower Hubbard subband. Con- 
sidering the arguments presented in HQ] we wm take into 
account the Coulomb interaction between the electrons lo- 
cated on the neighboring sites. The obtained system of the 
Hubbard fermions in the atomic representation will be de- 
scribed by the Hamiltonian of the t — U-model: 

h = E £ ° x r + E tf^TXrn + \ E w+a- 
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Here the first term reflects an ensemble of noninteracting 
electrons in the Wannier representation. The occurrence a 
fermion with spin projection a at site / increases the en- 
ergy of the system by value en, X^ q = \f,p){f,q\ are 



the Hubbard operators |17L 1 1 8ll describing the transition 
from the one-site state \q) to the state \p). The second term 
corresponds to the kinetic energy of the Hubbard fermions 
where matrix element t f m determines the intensity of elec- 
tron hopping from site / to site m. The last term of the 
Hamiltonian takes into account the Coulomb interaction 
between the electrons located on neighboring sites / and 
f + 5 with intensity V. The operator of the number of 
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electrons on site / is fij = XJ C T . 

a 

Below we limit our consideration to the case when the 
number of holes h = (1/N) ^/(^/°) m tne system is 
small; i. e., the condition h = 1 — n <S 1 is satisfied. 
In this regime it is reasonable to extract in explicit form 
the mean-field effects caused by the intersite interactions. 
Using the condition of completeness of the diagonal X- 
operators in the reduced Hilbert space Xj° +Xj^ +X^ = 
1 we express Hamiltonian £Q as 

H = -E + £( e - Wb)Xy + £ tfmX^X ^ 



So 
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+lE V (X°f -h)(X° f + s-h), 
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where E = 2NV(1 — h) 2 in the mean-field approxima- 
tion determines the energy of the Coulomb interaction of 
the system containing h holes per site. At h = value E 
equals to the exact value of the energy of the ground state 
of the system (with disregard for e ), since in this case hop- 
pings make no contribution. The renormalized value of a 
one-electron level e = e + 4V is determined by the fact 
that when there is one electron on each neighboring site 
the excitation energy increases by 4V. The shift e — 4hV 
is related to the decrease in the Coulomb repulsion energy 
when the average number of holes in the system is nonzero. 
Note that such mean-field renormalizations of the one-site 
energies of electrons were used previously, for example, 
in the Falicov-Kimball model during investigation of the 
transitions with the change in valence. Extraction of the 
obvious mean-field effects is needed for representation of 
the intersite interaction in the form containing the corre- 
lation effects only. One can see that the last term of the 
Hamiltonian (UJ will contribute only at the presence of no- 
ticeable fluctuations of the occupation numbers, i. e., when 
the SICs are relatively large. 

The method used below for the description of the strong 
intersite interactions (SII) consists in generalization of the 
Hubbard idea 0]] for consideration of the intersite interac- 
tions lllql . It follows from the exact equation of motion for 
operator 

ijX°°= (e - 4Vh)X? - V J2 XTW+s ~ h) 
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that at large V one should correctly consider the contribu- 
tions related to the second term in the right part of equation 
(01. One of the ways of solving such a problem is to ex- 
tend the basis set of operators by means of inclusion of the 
operators 



in which uncoupling is forbidden. This operator describes 
the correlated process of annihilation of an electron on site 
/, since the result of the action of operator X^" 7 depends on 
the nearest configuration neighborhood of site /. Recall for 
comparison that in the fundamental Hubbard's work |l|] the 
basis was expanded by means of addition of a set of one- 
site operators af a hj 9 . 

Inclusion of new basis operators requires the equations 
of motion: 

i^-fr = (e - 2Vh - V)& f - AMI - hWX / 
dt 1 1 1 

+ Y t fm ((x™ + x?)x% + X?X%) (X™ 5 - h) 

mS 
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-V ]T Xf{Xf +5 - h){Xf +5i - h). (5) 

<W<5i) 

One can see that among the operators containing large pa- 
rameter V a three-site operator has occurred, which reflects 
the correlation effects related to the presence of two holes 
in the first coordination sphere of site /. If hole concentra- 
tion in the system is low, then the contribution of this term 
can be ignored. Calculation of the spectrum with regard of 
this three-center operator confirmed validity this approxi- 
mation. 

After the basis set of operators has been specified, the 
set of equations of motion is closed using the Zwanzig- 
Mori projection technique Il9l l2p|| . In the main approx- 
imation one can neglect the kinetic correlators occurring 
after projection and limit the consideration to spatially ho- 
mogeneous solutions. Then, for the Fourier transforms of 
the Green functions we obtain the closed system of two 
equations: 

(u - £k )«X kCT |xL» = l ~Y^- 7 k «*w|*L»> (6) 

(w - &{{^\xD) = -4Mi - % k «x kCT |xL»- 

Here we made the following notation: 

£u = (s - AVh) + ( i±* V k , 7k = y_^S, 



& = (e - 4V7i) - V(l - 2h) + 
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*k = ^lk + At' cos k x cos k y + 2t"(cos 2k x + cos 2k y ), 
iik = 2t(cos k x + cos k y ). (7) 

The use of the spectral theorem allows obtaining the equa- 
tion which determines the dependence of the chemical po- 
tential on doping: 
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where f(x) = (exp(^^) + is the Fermi-Dirac func- 
tion, fi is the chemical potential of the system, and the two- 
band Fermi excitation spectrum is determined as 



= o T 



(9) 



+ 4h(i-/o 7 2. 



FiglU (on the right) demonstrates a band picture of the 
energy spectrum of the t — V-model obtained by solving 
equations ©. The values of hopping parameters were cho- 
sen so that evolution of the Fermi contour upon hole dop- 
ing would qualitatively correspond to that observed experi- 
mentally. The left part of the figure shows the energy spec- 
trum of the same model calculated with disregard of the 
SICs. Note that here the mean-field contribution of the SIIs 
is taken into account by the above-mentioned renormaliza- 
tions. Comparison of the two presented spectra shows that 
the correct account for the SICs yields the qualitative dif- 
ference, specifically, the occurrence of an additional band 
(BFS) in the band structure of the t — V— model. 

It is seen from (O that the resulting energy spectrum 
forms by hybridization of states from an ordinary Hubbard 
band with energies e k and the states induced by the fluc- 
tuations of configuration neighborhood with energies £ k . 
Note that unlike e k the dependence of function £ k on the 
quasimomentum is determined only by the integral of hop- 
ping between the nearest neighbors i lk . At V 2> ji k | the 
values of this function form an energy band located lower 
by V. The intensity of hybridization is determined by the 
value proportional to the one-site mean-square fluctuation 
of the occupation numbers 



(An) 2 = {(X° f ° - h)(X° f ° - h)) = h(l - h). (10) 

As a result the spectral intensity of the BFS in the regime 
7 » |ik| acquires the simple form: 



A-(k,w) 



4(l + /i) (An) 



1 + 16(An) 2 + V 1 + 16(An) 



(11) 

It follows from this formula that with an increase in h the 
spectral weight of the BFS rapidly grows and at h = 0.2 
the relative contribution of the BFS reaches 30%. 

The occurrence of the BFS qualitatively changes the 
electron density of states of the Hubbard fermions. The up- 
per graph in Fig|2]shows the density of states of the t — V- 
model with disregard of the SICs. It is seen that ignoring 
these correlations leads to trivial displacement of the band 
position. 

The energy dependence of the sum electron density of 
states of the t — V-model for the same set of parameters, 
but with regard of the SICs, is shown by a bold solid line 
in the middle graph of Figj2] Comparison with the upper 
graph evidences that the account for the SICs yields notice- 
able qualitative changes in the density of states: the latter 
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FIG. 1. Band picture of the t — V-model without account for 
the SICs (left) and with account for the SICs (right) for the set of 
parameters t' = —0.1, t" — —0.5, V — 2.5 (in terms of \t\) and 
holes concentration h — 0.15. Dashed lines show the chemical 
potential. 




FIG. 2. Density of states of the t — V-model calculated for h = 
0.25 without account for the SICs V = 1 (top), with account 
for the SICs at V = 1 (middle) an at V = 2.5 (bottom). Here 
t = —0.1, t" = —0.5 (all in terms of \t\). Dashed lines show the 
chemical potential. 



acquires a three -peak structure instead of the two-peak one. 
The occurrence of the additional peak corresponds to the 
formation of the BFS. For clarity, BFS density g~(E) is 
reflected in the graph by the line which bounds the shaded 
area. The BFS fraction is 36% of the whole number of 
states of the system, whereas the fraction of the ground 
band with density of states g + (E) is 64%. Summation of 
densities g + (E) and g~~(E) gives total density g(E) with 
the three peaks. 

More substantial changes in the energy structure of the 
t — I/-model take place at high values of the intersite 
Coulomb interaction. It is seen from the lower graph in 
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FIG. 3. Fermi surfaces calculated for the set of parameters t' = —0.1, t" = —0.5, V = 2.5 ( in terns of \t\) and different h without 
account for the SICs (dashed lines) and with account for the SICs (solid lines). 



FigJZl which presents the same dependences as in the two 
upper graphs but calculated for V = 2.5 (the rest param- 
eters are not changed). In this case, the SICs lead to split- 
ting the BFS off and the formation of a gap in the energy 
spectrum. The graphs presented in Fig|2] are calculated at 
h = 0.25. In the case of h = the SIC contribution 
becomes zero and the density of states of the system will 
be such as shown in the upper graph in Fig|2] This sug- 
gests the presence of a qualitatively new effect related to 
the account for the SIC: upon doping not only the chemical 
potential shifts but the density of states rearranges. 

The growth of the BFS spectral intensity with an increase 
in doping level leads to renormalization of the dependence 
of the chemical potential on hole concentration in the sys- 
tem. As a result in the area of optimal doping h ~ 0.2 
the square bounded by the Fermi contour noticeably grows. 
The Fermi contour calculated at h = 0.25 with disregard 
of the SICs is shown by dotted lines in Fig(3] If the cor- 
relations are taken into account then the Fermi contour no- 
ticeably grows and acquires the form shown by solid lines. 
In this case the change in the square bounded by the Fermi 
contour reaches 16% (right graph). This increase is im- 
portant for description of the Lifshitz quantum phase tran- 
sitions occurring upon doping rf2~ill and for interpretation 
of the experimental data on measuring magnetic oscilla- 
tions in the de Haas-van Alphen effect. Recently, such 
measurements have been performed by many researchers 
due to substantial improvement of quality of materials and 
novel techniques with the use of the strong magnetic fields. 

Note in summary that the hybridization character of the 
energy spectrum with regard of the SICs is related to the 
fact that electron hoppings between the nearest neighbors 
lead to the transitions between the one- site energy levels 
which differ by V. Therefore, the intensity of such hy- 
bridization processes is proportional to both hopping pa- 
rameter t and the above-mentioned mean-square fluctua- 
tion of the occupation numbers. Note also that the consid- 
ered modification of the energy structure due to the SICs is 
general and not limited merely by the t — T/-model. One 
should expect that the effects discovered in this study will 
be especially important for the systems with variable va- 
lence. 
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